Proton Production in d+Au Collisions and the Cronin Effect 
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Proton production in the intermediate px region in d+Au collisions is studied in the parton recom- 
bination model. The recombination of soft and shower partons is shown to be important in central 
collisions, but negligible in peripheral collisions. It is found that the large nuclear modification 
factor for proton production can be well reproduced by a calculation of the 3-quark recombination 
process. 
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In a previous paper Q we have shown that the Cronin 
effect |2| on pion production in d+Au collisions can be 
understood in terms of the recombination of the soft and 
shower partons without px broadening by multiple scat- 
tering in the initial state. In this short note we extend 
the consideration to proton production and show that 
the same effect can similarly be interpreted. 

The Cronin effect can best be displayed by the nuclear 
modification factor Rcp(pt) that is the ratio of central- 
to-peripheral yields appropriately scaled by the average 
number of binary collisions (iVcoii)- It is found in the 
PHENIX experiment that R^p for proton reaches a value 
roughly 2 for 2 <Px < 3 GeV/c, even higher than that 
for pion at ~ 1.4 Such a behavior of the enhancement 
effect is hard to interpret, if hadrons produced at inter- 
mediate pt are the consequences of fragmentation of hard 
partons produced at higher px- Indeed, since there is no 
energy loss in d+Au collision, one would expect Rcp ~ 1 
for both pion and proton on the grounds that fragmenta- 
tion outside the cold medium should be independent of 
the impact parameter. Thus the observed R^'p strongly 
suggests the dependence of the hadronization mechanism 
on the medium. In our view recombination is that mecha- 
nism, which, on the one hand, provides a way to describe 
fragmentation in terms of shower partons , and, on the 
other, can take into account the coalescence of soft and 
shower partons to form hadrons in the intermediate pt 
range [5j- 

For the proton spectrum at Iow^t, one should be care- 
ful about the mass effect. The low-pr region is, however, 
not the main part of our work where the model has any 
predictive power. As in P, we fit the data in that 
region and use our model to predict the hadronic spectra 
at intermediate and high px- We have found that for 
the purpose of data fitting at low px our ID formula- 
tion of the recombination process is quite adequate when 
only the kinematical variables are suitably modified to 
account for the proton mass. Thus we start with the 
invariant inclusive distribution for proton formation at 
midrapidity in the recombination model 0, ||| 
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where all momentum variables pi and p are transverse 
momenta, andp denotes the energy of the proton. Since 



the parton masses are set to zero, we continue to use 
Pi for their energies. F(pi,p2,P3) is the joint distribu- 
tion of it, it, and d quarks at pi,P2 and P3, respectively. 
Rp(Pi>P2,P3,p) is the recombination function for a pro- 
ton with momentum p [fj 0] 
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where yi — Pi/p, a — 1.75, 7 = 1.05, and 

g= [6B(a + l,a + 7 + 2)B(a + l, 7 + l)r 1 

B(a, b) being the beta function. As always in the recom- 
bination model, the main issue is about the distribution 
of the quarks that recombine. Here, it is F(pi,p2,P3). 

Following the same notation used in p| for Au+Au 
collisions, we write schematically 



F = TTT + TTS + TSS + SSS, 



(4) 



where all the shower partons iS are from one hard par- 
ton jet. Shower partons from different jets are ignored 
here for RHIC energies. T denotes thermal parton, even 
though in rf+Au collisions the notion of thermal equilib- 
rium may not be justified. To preserve the same notation 
as in , we continue to use T to signify the soft partons 
that are not associated with the shower components of a 
hard parton and are loosely referred to as thermal par- 
tons when convenient. The SSS term in Eq. (@J), when 
convoluted with R p in Eq. QJ, gives rise to what is usu- 
ally regarded as the fragmentation of a hard parton into 
a proton jj| . The TTT term comes entirely from the soft 
partons, while TTS and TSS accounts for the interplay 
between the thermal (or soft) and shower partons. 

As in P, Q , we can calculate the distributions of the 
shower partons S from the QCD processes of producing 
hard partons and their induced shower partons. We can- 
not calculate the thermal component T, which is deduced 
from fitting the low-px data. Thus what we can calculate 
that is new is only the effect of the TTS + TSS terms at 
intermediate px, knowing that the <S<SiS term dominates 
at very high px ■ 

For T we use the same parametrization as before P, || , 
and write 



T(pi) =pi 



dNj 



dpi 



Cpi exp(-pi/T), 
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where T should be regarded as just an inverse slope. The 
thermal contribution to the proton spectrum arising from 
TTT recombination is then 

dtC = C^p^ _ p/T B(a + 2,7 + 2)B(a + 2, a + 7 + 4) 
pdp 6 p° e B(a + 1,7+ l)B(a+ l,a + 7 + 2)'^ ' 

which differs from a similar formula in |5J] by only the 
presence of p°, instead of p. For the other three terms 
in Eq. (J3J that involve S, the contributions to the pro- 
ton spectrum are the same as those in jjj except that 
each equation should be multiplied by the factor p/p° on 
the RHS and the factor £ should be omitted. The lat- 
ter refers to the mean energy loss in Au+Au collisions, 
and should be absent in d+Au collisions. Also the hard 
parton distributions fi(k) in J| should be changed to the 
corresponding ones for d+Au collisions, as given in 0. 
What is to be emphasized is that there are no free pa- 
rameters to adjust in those terms. The shower parton 
distributions are the main input, and they have previ- 
ously been determined in 

We must now determine C and T by fitting the low- 
Pt data using Eq. JBJ for both central and peripheral 
collisions. The data available are from PHENIX, given 
as figures online We fit them in the region 0.5 < 
Pt < 1.5 GeV/c and the results are shown by the light 
solid lines in Fig. 1 (a) and (b) for 0-20% and 60-90% 
centralities, respectively. The values determined are C = 
11.5 (8.0) (GeV/c)" 1 and T = 0.24 (0.21) GeV/c for 0- 
20% (60-90%) centrality. The contributions from TTS + 
TSS and SSS components are determined without free 
parameters, and are shown by the dashed and dash-dot 
lines in the same figures. The sum of all four components 
are shown by the thick solid lines. Evidently, they agree 
with the data very well. 

We note that at 0-20% centrality the thermal-shower 
(TTS + TSS) contribution crosses over the fragmenta- 
tion (SSS) component at px ~ 2.5 GeV/c, roughly the 
same as in the case of pion production in d+Au colli- 
sions p|. However, the thermal TTT contribution is 
roughly the same as each of the (TTS + TSS) and SSS 
contributions at the cross-over point, whereas for pion 
production the thermal TT contribution is much lower 
than TS and SS at the same point. Thus the thermal 
contribution to proton formation dominates over a wider 
range of px than that for pion. That is because of the C 3 
dependence in Eq. (JSJ), and is consistent with the find- 
ings in Refs. [iiUJl, where the recombination of thermal 
partons can account for the large p/ir ratio up to pt ~ 
3-4 GeV/c Au+Au collisions. The same cannot be said 
about 60-90% centrality in d+Au collisions. The cross- 
over between TTS + TSS and SSS occurs at pt ~ 1 
GeV/c, where the distributions are far lower than the 
thermal contribution. Throughout allpx in Fig. 1 (b) the 
(TTS + TSS) component is negligible, a situation that is 
similar to what we have found separately for either p or tt 
production in pp collisions. The reason is the substantial 
reduction of the soft parton environment in peripheral 
d+Au or pp collisions so that the (TTS + TSS) compo- 
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FIG. 1: (a) Proton transverse momentum distribution in 
d+Au collisions at 0-20% centrality. The preliminary data 
are from |8(. (b) Same as in (a), but for 60-90% centrality. 



nent does not develop any strength relative to the SSS 
component. We should, however, note parenthetically 
that the low density of soft partons that are produced 
in pp collisions is sufficient to make the jet structure dif- 
ferent from that in e + e~ collisions, which have no soft 
partons of the same origin at all. 

The fact that our calculated results agree well with the 
data for both 0-20% and 60-90% centralities implies that 
the ratio BF cp (pt), defined by 



R p CP (pi 



(^Coii) 60 _ 90% dN p /d PT (0 - 20%) 
(^Coii} _ 20% dN p /dp T (m - 90%) 



(7) 



must also agree with the data. That agreement is shown 
explicitly in Fig. 2 by the solid line, where the data are 
from Ref. Q ■ The calculated curve approaches 1 at high 
Pt, where the yields are dominated by the fragmenta- 
tion of hard partons (SSS), which is independent of the 
soft partons. The dashed line in Fig. 2 shows the contri- 
bution to R P CP from the thermal components only. The 
small difference in T for the two centralities results in 
an exponential growth in p? as can be seen directly from 
Eq. jnj). Thus the effect of thermal-shower recombination 
is the damping of that exponential increase in R p -, p (px), 
as shown by the solid line. Because of the ineffectiveness 
of the thermal-shower contribution at 60-90% centrality, 
that damping does not take place until px reaches near 2 
GeV/c. By then R P CP is already greater than 1.7, which 
is higher than Rq P . Thus the origin of R^ P being greater 
than Rq P is mainly in the TTT recombination for pro- 
ton being more sensitive to centrality than TT for pion. 
The role of shower partons is limited in that comparison. 
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It should be pointed out that our fitting procedure in 
the determination of C and T has not ignored R^ P as 
an outcome. Since Fig. 1 has logarithmic vertical scale, 
those data points can determine C and T only within 
narrow ranges. The data in Fig. 2 involve their ratios 
and are in linear scale. Thus the parameters can be more 
accurately determined by including Fig. 2 in the fit. Since 
the distributions of the soft partons at low px are not 
generated from first-principle calculations, we have taken 
the liberty to use all the data available to make the best 
determination of them. Our predictions are only for the 
behavior at px above the region that is dominated by the 
soft partons. 

The values of C and T obtained here should be com- 
pared with those determined from the pion spectrum Q . 
They are = 12 (5.65) (GeV/c)- 1 and T« = 0.21 (0.21) 
GeV/c for 0-20% (60-90%) centrality, where the sub- 
scripts 7r are added for distinction. Whereas there is 
no essential dependence of some of those parameters 
on whether the formed hadrons are pions or protons, 
i.e., C p ~ C n at 0-20% and T p = T n at 60-90%, there 
exist some significant differences in others: C v = 8.0, 
C„ = 5.65 at 60-90% and T p = 0.24, T„ = 0.21 at 0-20%. 
The species dependence of those parameters reflects the 
general properties of the spectra, especially at px < 1 
GeV/c. In [H it is shown in Au+Au collisions that 
the low-pT spectra can be fitted by exponential behav- 
ior in m,T with the inverse slope increasing linearly with 
hadronic mass at a rate that increases with centrality. 
It strongly suggests hydrodynamical expansion radially, 
which is well known to exist. It means that a significant 
portion of the hydrodynamical fluid is hadronic at very 
low pt- Only the portion that remains partonic at pt 
around 1 GeV/c and above are available for recombina- 
tion with the shower partons. How much of this scenario 
is valid (and in what quantitative way) for d+Au colli- 
sions is not known. Since our hadronization model does 
not treat the hydrodynamical expansion phase, we can 
only take the thermal hadrons in the 0.5 < px < 1-5 
GeV/c range as observed, but not lower, to determine 
our parameters for the thermal partons T. The species 
dependence of some of the C and T parameters is clearly 
a result of our inadequacy in subtracting out the already- 
formed hadrons from the medium at low px, since the 
thermal partons T should have no knowledge of what 
hadrons they are to form. Given our inability to treat 
very low px physics, we can only regard what we have 
done as a determination of the intermediate px behavior 
in the separate cases of specific hadrons without a way 
to enlighten the problem of overall species dependence at 
very low px in a broader scheme. 

The important properties of hadron production in 
d+Au collisions that we have learned from this study is 



that the protons are formed by recombination at all px 
and that the underlying partons that give rise to their for- 
mation change smoothly from the soft component to the 
semi-hard shower partons that are created by the hard 
partons. The recombination formalism allows us to calcu- 




FIG. 2: Nuclear modification factor Rq P for proton produc- 
tion in d+Au collisions. The data are for 0-20% to 60-90% 
centralities The solid line is the result of our calcula- 
tion when all contributions are taken into account, while the 
dashed line gives the ratio when only the thermal contribu- 
tions are included. 



late the px distribution in the intermediate and higher px 
regions with good agreement with the data. The contri- 
bution from the recombination of soft and shower partons 
cannot be interpreted as a modification of the fragmenta- 
tion function, since the hard partons in d+Au collisions 
are not significantly affected by the cold medium that 
they traverse. In this treatment the Cronin effect of the 
proton spectrum is not induced by transverse broaden- 
ing due to initial scatterings, but is caused mainly by the 
centrality dependence of the soft partons that recombine. 
Rqp is higher than Rq P because the number of such re- 
combining quarks is 3 instead of 2. Our approach cannot 
be viewed as being totally successful until the data for 
px > 3 GeV/c turn out to support our predictions in the 
higher px region. 
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